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1. Introduction 

Precise determination of the CKM matrix is an important step to establish the validity range of 
the Standard Model, thereby setting the basis for exploring physics at smaller space-time scales. In 
this respect the CKM matrix elements in the second row, especially \V c d\ and |V CS |, are still poorly 
determined with errors in the 10% range [JTJ] _ This situation contrasts sharply with those of the first 
row where \V u d \ and |V Mi | are known at sub percent level. 

Lattice QCD provides a variety of means for precise determination of the second row of the 
CKM matrix. With \V CS \, for example, the leptonic decay width of the D s meson T(D S — > Iv) is 
given by 

T{D s ^lv) = ^m}m D jl s (^-^\ \V CS \ 2 . (1.1) 

so that a lattice determination of the decay constant fjj t combined with the experimental value of 
T(D S —> Iv) will allow us to extract \V CS \. The other matrix element |V C ^| can be obtained in the 
same way. 

An intriguing topic in the charm quark sector is a disagreement of the value of the D s meson 
decay constant from lattice QCD and experiment. In 2007 HPQCD and UKQCD Collaboration 
reported a calculation of this quantity using the HISQ action for heavy quark on the 2+1 dynamical 
flavor MILC configurations [0]. Their value fn s = 241 ± 3 MeV exhibited significant deviation 
from those of the CLEO and Belle experiments at the time. The latest CLEO value fo s = 
259.5 ±7.3 MeV, though smaller than 2007, is still two standard deviations away. 

FNAL Collaboration also calculated the decay constant on the MILC configurations using the 
clover heavy quark formalism. Their original value fp s = 249(3)(16) MeV [Q] was also smaller 
than experiment, but an update this year [^] including a 2% change in the scale setting of the MILC 
ensemble moved the value up to fo s = 260(10) MeV. 

In this paper we present a status report of our work on the charmed mesons using the relativistic 
heavy quark formalism on the 2+1 dynamical flavor PACS-CS configurations. 



2. Setup 

The charm quark system is simulated with a relativistic heavy quark action on the 2+1 flavor 
lattice QCD configurations. The configurations are generated by the PACS-CS Collaboration [||] 
using the nonperturbatively 0(a)-improved Wilson quark action with = 1.715 ^ and the 
Iwasaki gauge action. The lattice size is 32 3 x 64 whose spatial extent is L = 3 fm with the lattice 
spacing of a = 0.09 fm. The dynamical up-down quark masses range from 8 MeV down to 3 
MeV, which is close to the the physical value. We utilize the same quark action for the light quark 
propagators in our calculation. 

For the heavy quarks, we employ a relativistic heavy quark action proposed in Ref. ^J. The 
cutoff errors are reduced from 0{{niQa) n ) to 0(/(mga)(aAgcD) 2 ) where fipiQa) is an analytic 
function around niQa = 0. We use one-loop perturbative values for parameters in the heavy quark 
action [^]. In addition, nonperturbative contributions at the massless limit is included for the clover 
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Table 1: Simulation parameters. Quark masses are perturbatively renormalized in the MS scheme. The 
renormalization scale is ji = I /a for each simulation point and ji =2 GeV for the physical point. 
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Figure 1: m^f quark mass dependence of mo(left panel) and mu s (right panel). 



coefficients, and the parameter v is adjusted nonperturbatively to reproduce the relativistic disper- 
sion relation of the charmonium. Our valence heavy quark masses are always tuned to the physical 
charm quark mass point, which is determined by the spin-averaged IS state of the charmonium, 
M(15) = {M^ c +3M 7/ y)/4 = 3.0677(3) GeV. Our results for charmonium spectrum have been 



reported in Ref. [ 1C] 



Table | summarizes our simulation parameters and the statistics of the configuration sets we 
have used for the heavy quark measurements. The number of the source points is quadrupled to 
suppress statistical fluctuations. 

3. Results for D and D s meson masses and decay constants 

We present our results for D and D s meson masses and their decay constants. Since the physi- 
cal charm quark mass has already been fixed with the heavy-heavy spectrum, all heavy-light meson 
quantities can be predicted. 

Figure [I] compares our results for the D and D s meson masses with the experimental values [||]. 
Since our sea quark masses are close to the physical values, we extrapolate our results with a linear 
function of the up-down and the strange quark masses to the physical point, 

M = A + B(m ud - m^ ys ) + C(m s - mf? 8 ). (3. 1) 
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Figure 2: quark mass dependence of /b(left panel) and fo s (right panel). 
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Figure 3: m^? quark mass dependence of fo s over fy. 



Our results at the physical point (red bursts) reproduce the experimental spectrum well. 
We also calculate the decay constants of heavy-light pseudoscalar mesons, 

<0\A*? p \PS(p)> = ifpsPi l , (3-2) 
A™ p = y/2%y/2^Z All {q(x)Ypy$Q(x) (3.3) 

-c+d+ {q{x) W Q{x)) - cjd- {q{x) W Q{x)) } . (3.4) 

For the renormalization factor and the improvement coefficients of the axial current, we employ 
one-loop perturbative values [11]. Furthermore, the nonperturbative contribution at the massless 
limit is incorporated to the improvement coefficient c\ by writing 

ct = (cl(m Q a)-cl(0)r + cT (3.5) 

with cf = -0.03876106 0. 

Figure |2| shows the decay constants, where we also plot other 2+1 flavor lattice QCD results 
from HPQCD and UKQCD Collaboration ^ and FNAL group |5[], as well as the experimental val- 
ues [p], [3"| |I3[], for comparison. Our value for fo s agrees with the experimental determinations, while 
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Figure 4: quark mass dependence of hyperfine splittings for charm-Mii(left panel) and and charm- 
strange mesons (right panel). 



that for fo is somewhat larger. Comparing the three sets of lattice determinations, we observe, both 
for fo and fo s , an agreement between our values and those of the FNAL group, while there seems 
to be a clear discrepancy between our values and that by the HPQCD and UKQCD Collaboration. 
Possible reasons are our use of perturbative renormalization factors and the necessity of taking the 
continuum extrapolation. 

In Fig. H we plot the ratios of fu s to fo in which uncertainties coming from the perturbative 
renormalization factors cancel out, and perhaps that of the lattice cutoff to some extent. It is 
interesting that all three lattice results are mutually consistent with small errors of a few percent, 
but that they seem to lie at the lower edge of the experimental value. 

4. Results of hyperfine splittings for charm- ud and charm-strange mesons 

We evaluate the hyperfine splittings for charm-i«i and charm-strange mesons. Masses of D* 
and D* are obtained from their two-point functions. D* and D* decays are prohibited energetically 
on our lattice. 

Figure [| represents our results for the hyperfine splittings with the experimental values [jj]]. 
For D*-D mass difference, our result at the physical point (red bursts) reproduces the experimental 
spectrum. On the other hand, D*-D s mass difference shows 10% deviation from the experimental 
value. A possible origin of this discrepancy is the scaling violation. 

5. Estimating the CKM matrix elements 

We attempt to estimate the CKM matrix elements from our data of D and D s meson masses 
and decay constants combined with experimental values of leptonic decay widths. Using the CLEO 
value of T(D S -► Zv) [g], we find 

\V cs \{lattice) = 0.98(2)(3) + 0(g 2 a), (5.1) 
\V CS \(PDG) = 1.04(6) § (5.2) 
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The first error is statistical, and the second is experimental. For comparison, the PDG value is 
listed. Our preliminary result is consistent with the PDG value. We note that our estimate of |V CS | 
still has a discretization error of 0(g 2 a) due to the use of perturbative renormalization factor. We 
must take the continuum extrapolation to remove this source of uncertainty. 

We can estimate \V c d\ from the D meson mass and decay constant and CLEO value of F(D — > 
IV) ©: 

\V cd \(lattice) = 0. 201 (2) (9) + 0(g 2 a), (5.3) 
\V cd \(PDG) = 0.230(11) (5.4) 

Our estimate is smaller than the PDG value by about 10%. For completeness we also record the 
ratio |V fS |/|V CY /| for which 0(g 2 a) errors drop out, and is replaced with scaling violation of 0(a ): 

^(lattice) = 4.72(13)(26) + <9(a 2 ), (5.5) 

(PDG) = 4.52(34) 0. (5.6) 



\Vcd\ 
\V cd 



6. Conclusion 



We measured charm-w<i and charm-strange meson masses and decay constants with the rela- 
tivistic heavy quark action on the 2+1 dynamical flavor PACS-CS configurations at cT x = 2.2 GeV. 
Since our sea quark masses are close to the physical point, sea quark mass corrections in our results 
are small and only short extrapolations to the physical point are needed. 

We found our results are consistent with the experimental values in two standard deviations. 
We do not observe discrepancy from experiments in fu s . 

Combining our data of masses and decay constants with experimental values of leptonic decay 
widths, CKM matrix elements are evaluated. Our estimates still have discretization errors, which 
is expected to be 1% by the order counting. Continuum extrapolation is needed to achieve a few 
percent accuracy. 
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